Unconventional metal-insulator transition in two dimensions 
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We show, by using a correlated Jastrow wave function and a mapping onto a classical model, that 
the two-dimensional Mott transition in a simple half-filled one-band model can be unconventional 
and very similar to the binding-unbinding Kosterlitz-Thouless transition of vortices and anti- vortices, 
here identified by empty and doubly occupied sites. Within this framework, electrons strongly 
interact with collective plasmon excitations that induce anomalous critical properties on both sides 
of the transition. In particular, the insulating phase is characterized by a singular power-law behavior 
in the photoemission spectrum, that can be continuously connected to the fully-projected insulating 
state, relevant to strongly correlated low-energy models. 
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I. INTRODUCTION 

The metal-insulator transition (MIT) driven by elec- 
tron interaction, the so-called Mott transition^ is one 
of the most challenging issues in modern solid state 
physics, especially because of its possible connections 
with other phenomena, like high-temperature supercon- 
ductivity. The prototypical model to study the MIT is 
the one-band Hubbard model, where the transition shows 
up if the ratio between the bandwidth W and the on-site 
Coulomb repulsion U is varied^ In a pioneering work, 
Brinkman and Rice 3 argued that the MIT could be con- 
tinuous, however, their approach led to a Mott insulator 
where charge fluctuations are completely frozen. Recent 
developments, based on dynamical mean-field theory cal- 
culations^ are able to give a more realistic description of 
the MIT where, at zero temperature, a wide gap in the 
excitation spectrum opens immediately in the insulator. 
However, within this approach, the role of the dimen- 
sionality is not taken into account, whereas its relevance 
comes out from recent experiments on Si MOSFET's 5 
and organic materials^ Moreover, the situation can be 
radically different whenever the Hamiltonian contains a 
true long-range Coulomb interaction, as pointed out in 
the original Mott argument. 1 

In this work, we focus our attention on a correlated 
wave function (WF), which is expected to closely de- 
scribe the MIT in two dimensional (2D) systems with 
long-range Coulomb interactions. The important fact 
is that, in the strong-coupling regime, we can adiabat- 
ically connect our WF with the fully-projected one, usu- 
ally considered to describe systems in the limit of in- 
finite Coulomb repulsion, e.g., the so-called resonating 
valence-bond stateii§i2ii£ Moreover, the metallic phase 
has no quasiparticles defined, showing low-energy prop- 
erties similar to the one-dimensional Luttinger liquid. 

The paper is organized as follow: In Sec. [H] we show 
the mapping between the quantum wave function and a 
classical model at finite temperature and in Sec. II I II we 
show our numerical results and we draw the conclusions. 



II. CLASSICAL MAPPING 

Let us discuss how to construct a WF for correlated 
insulators. In general, starting from the ground state 
|\&o) of a system with N electrons, with energy Eq, it 
is possible to construct simple variational states for the 
lowest-energy excitations. For instance, in analogy with 
the Feynman's construction for the liquid Helium^i the 
plasmon excitation with momentum q is given by: 
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where n q is the Fourier transform of the local electron 
density. Its variational energy is 

(-k)q 2 



Eq = E 



2N„ 



(2) 



where (k) is the ground-state kinetic energy per particle 
and 



N q = 



{^a\n q n^ q \^ ) 
(*o|*o) 



(3) 



is its static charge structure factor. When applied to an 
insulator, with gapped charge excitations, this implies 
that, for small momenta, N q ~ q 2 . This argument has 
very general consequences on the form of \^o), that do 
not depend on the particular microscopic model. To this 
purpose, let us denote an electronic configuration by the 
positions {x} of the particles. For all the operators 9 that 
depend only on such positions, e.g., the structure factor 
itself, the quantum average 



(0) = 



(*o|*o) 
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can be written in terms of the classical distribution 



|*o(x)| 2 = |( S |*o}| 2 /E^I(^l*o>| 2 ,as 

(e) = J2(xm\Mx)\ 2 



(5) 
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Since 1^0(^)1 2 is a positive quantity, we can define an 
appropriate correspondence between the WF and an ef- 
fective potential V(x): 



2 = e -V(, 



(6) 



The size consistency of the WF implies that the poten- 
tial V(x) is extensive, namely of order N for typical con- 
figurations. In the limit of strong Coulomb interactions, 
there are small charge fluctuations and, therefore, we can 
safely assume that only the two-body term is relevant and 
all multi-particle interactions are negligible. This leads 
to the quadratic potential 



q^O 



(7) 



being n q {x) the Fourier transform of the local density of 
the configuration \x). To obtain the expected behavior 



N q = S ^n q (x)n_ q (x)t 



-V(x) 



the effective potential must diverge as 

— T-j—^ + less singular terms. 
T e} j q 1 



„eff 



(8) 
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Here T e ff can be considered as the effective temperature 
of classical charges interacting through a potential n/q 2 . 
Within this choice of v q ^ , N q ~ q 2 is generally valid 
and can be understood by considering n q as a complex 
continuous variable, so that the classical average oln q n- q 
turns into a standard Gaussian integral, yielding 



N„ 



1 



T eff q 2 



,eff 



(10) 



It should be noted that the fully-projected wave function 
with no charge fluctuations, and therefore N q = for 
\q\ ^ 0, is recovered when T e ^ — > 0. 



III. RESULTS 

Let us now consider a general one-band fermionic sys- 
tem in 2D, in which every site of a square lattice can be 
either empty, singly occupied, being the electron either 
with spin up or down, or doubly occupied. A true Mott 
insulator, that does not break any lattice symmetry, can- 
not be represented by a simple WF containing a single 
determinant, and, at this stage, it is useful to define a 
state that is simple enough and yet is compatible with 
the predicted form of v*** in the insulating phase. A 
straightforward way to modify the effective potential de- 
termined by an uncorrelated determinant \T>) is obtained 
by taking into account an appropriate Jastrow factor J: 



\9) = J\V), 



(11) 
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FIG. 1: (Color online) Inverse of the dielectric function 1/e 
[see Eq. 1141 ] for the free-electron determinant. Left panel: 
1/e as a function of the effective temperature 1/(3 and for 
different sizes L of the cluster. The critical temperature of the 
classical Coulomb gas model y c CGM [ s marked with a dashed 
line for a comparison. Right panel: size scaling of 1/e for 
various (3. 



here \T>) is an electronic determinant that will be speci- 
fied in the following and J is a Jastrow term that depends 
upon the electronic density: 



J — exp 



q^0 



(12) 



where v q is the Jastrow potential, whose small-g behavior 
is given by: 



7T/3 



2 [2 — (cos q x + cos q y )} 
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(3 fixing its strength. 

At half filling, in onc-dimcnsional electronic systems, 
we founcUSiii that the singular Jastrow v q ~ ir(3/q 2 al- 
ways leads to an insulator, for any positive 0. In 2D the 
situation is different and a much more interesting sce- 
nario is obtained, with a phase transition as a function 
of the correlation strength (3. Indeed, given the behav- 
ior of the Jastrow v q ~ Tr(3/q 2 , the potential V(x) of 
Eq. 10 turns out to be the one of the classical Coulomb 
gas model (CGM). In this approach, particles with charge 
qi, corresponding to empty (qt — 1) and doubly = —1) 
occupied sites, interact through a Coulomb potential in a 
neutral background, represented by singly occupied sites 
(qi =0). In the half-filled case, there is an equal number 
of empty and doubly occupied sites, implying the charge 
neutrality of the CGM. The fugacity z of the CGM, that 
sets the average value of the charges, can be identified 
with the on-site Gutzwillcr term in the Jastrow potential, 
i.e., z = exp(— g), where g is the Gutzwiller parameter. 

The classical CGM in 2D is particularly interesting be- 
cause it shows a Kosterlitz-Thouless (KT) transition at a 
finite temperature T^ GM M> This transition is transpar- 
ent from the classical dielectric function: 
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FIG. 2: (Color online) The same as in Fig. Q for a gapless 
BCS state with A/t — 1.1 and d x i_ y i symmetry. 



where T e ^ is the temperature of the classical model. 
The charge structure factor is quadratic at small mo- 
menta, i.e., N q ~ aq 2 , for all temperatures, but the coeffi- 
cient a changes discontinuously at T GGM . Above T GGM , 
the CGM is in the plasma phase, i.e., a metallic phase 
with infinite dielectric function, perfect screening, and 
exponential correlation functions. On the other hand, 
below T GGM , the CGM is in the confined phase, with a 
finite dielectric constant. In this phase the charges are 
bound together forming dipoles, that, because of their 
residual interaction, induce power-law correlations. At 
the transition, the inverse of the dielectric function has a 
finite jump, changing from zero, in the plasma phase, to 
a finite value, in the confined phase. 

A similar mapping between a quantum state and a clas- 
sical model has been emphasized also in the context of 
the fractional quantum Hall effect: the Laughlin WF can 
be related to a classical system with particles interact- 
ing through a logarithmic potential. 15 However, in this 
case, all the particles have the same charge, forming a 
one-component plasma, and by varying the strength of 
the potential, there is a first-order transition between an 
incompressible fluid and a Wigner crystal^ The pecu- 
liarity of our approach is that, due to the mapping onto 
the two-component CGM, it is possible to connect con- 
tinuously the plasma phase to the insulating one. 

In the following, we show that, in analogy with the 
classical CGM, also in the case of fermionic systems at 
zero temperature, a KT-like transition is found by vary- 
ing the correlation strength that can be tuned by the Jas- 
trow strength [3. However, the existence of the fermionic 
part induces non-trivial properties for the two phases in- 
volved, that are not present in the classical problem. For 
example the uncorrelated part of the WF may contribute 
to the expression of the effective temperature T e ", as 
shown below. Whenever the square of the WF describes 
the plasma phase of the corresponding classical model, 
we can safely assume that the Gaussian fluctuations are 
exact for small g's and the classical temperature can be 
determined by imposing 1/e = in Eq. I|14fl . namely 
T e ff = 2ir lim 9 ^o N q /q 2 . In the language of quantum 
states, the Gaussian approximation leads to the well- 
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FIG. 3: (Color online) Equal-time density structure factor 
N q for the correlated wave function of Eq. (II IB (full squares) , 
compared to the same quantity calculated within the Gaus- 
sian approximation [indicated as GSA and given by Eq. 115H 1 , 
(full triangles) for [3 = 4/tt (upper panel) and [3 — 12 /it (lower 
panel, notice the different scale of the GSA data). 



known expression 



N q = 



N° 
q 



l + 2v q Nl 



(15) 



where N® is the charge structure factor of the uncor- 
related determinant |X>)iALi& The previous form of N q 
allows us to identify the effective temperature as: 



1 



T e ff 



(16) 



where ao = lim g ^ Q 2 /N® . 

In order to show the general validity of our approach, 
we consider the case of a free-electron determinant, ob- 
tained by occupying the lowest-energy states in the tight- 
binding model with dispersion Ek — — 2i(cos k x + cos k y ) 
and a gapless BCS state with a superconducting order pa- 
rameter A/; = A^osfcz — cos k y ). In these cases ao = 
and, therefore, the effective temperature in Eq. ll(j|) 
is determined only by the Jastrow coefficient, namely 
T e ff = 1/(3. In Fig. ^ we report the inverse of the 
dielectric function for the free-electron determinant and 
different sizes L of the system at half filling, i.e., N = L. 
In order to have closed-shell states for \V), we used 2D 
square lattices tilted by 45° (i.e., with L = 2l 2 and l x 
odd) and periodic boundary conditions. By increasing 
L, the curves show a steeper and steeper shape in the 
vicinity of the critical temperature T c . This result is fur- 
ther confirmed by the size scaling of 1/e, which clearly 
supports the existence of a finite jump in the thermody- 
namic limit: 1/e — > for T e ^ > T c , whereas 1/e — > const 
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FIG. 4: (Color online) Quasiparticle weight Zk 
(n/2,n/2) for the gapless BCS state with A/t = 
d x 2_ y 2 symmetry as a function of L and for different Jastrow 
strengths (3 (full circles). The case of the fully-projected wave 
function (empty circles) is also reported for A/t = 1.1 and 0.5. 



1.0 



0.8 



0.6 



0.4 



0.2 - 



X 



0.2 

D 



% 



5.0 10 U/t 15 ° 



• : L=98 
■: L=162 
: L=242 



5.0 



10.0 



Www 



u/t 



15.0 



FIG. 6: (Color online) Quasiparticle weight Zk at k = 
(7r/2, 7t/2) of the optimized paramagnetic wave function con- 
taining a Jastrow factor applied to the Fermi g func- 
tion of the interaction U /t in the Hubbard model, for three 
different sizes of the system. Inset: the number of double 
occupancies D as a function of U/t. 
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FIG. 5: (Color online) The behavior of 9 [the exponent of 
the quasiparticle weight, see Eq. 1191 1 as a function of j3 for 
BCS state (full circles) and the Fermi gas (FG) determinant 
(full squares). The value of the fully-projected states are also 
reported (arrows). 



for T e ^ < T c . Interestingly, T c depends slightly upon 
the choice of the uncorrelated determinant (see for com- 
parison Fig. |21 for the gapless BCS state) and is quite 
close to the CGM critical temperature T C CGM = 1/4. 
These results give an important and transparent insight 
into the strong-coupling limit described by the fully- 
projected WFjX that can be connected to our WF by 
letting j3 — > oo, i.e., T e ^ — > 0. Indeed, in the confined 
phase for T e *f < T c , the classical KT scaling equations 
of the CGM flow to fixed points with zero fugacity: this 
translates into the fact that the fully-projected state rep- 
resents the fixed-point of the correlated WFs describing 
the 2D Mott insulating phase. Therefore, in the con- 
fined regime, the ground-state properties are universal 
and represented by the ones of the fully-projected WF. 
In this respect, the total projection is not an unrealistic 
assumption and can accurately reproduce the low-energy 
physical properties of a strongly correlated system. On 
the other hand, for T e ^ > T c the classical KT scaling 
equations flow to strong coupling and are useful only close 
to the transition point. 

In Fig. |3 we show that in the plasma phase the Gaus- 



sian approximation, given by Eq. i|15|) . is very accurate, 
not only for small g's (where it is exact) but also for large 
momenta. In this case, the cusp singularity in N q for 
Q = (tt, 7r), related to the Friedel oscillations, is not re- 
moved, even though N q ~ q 2 at low momenta. Thus, for 
T e *f > T c , the WF (f[J describes a "Coulomb metal", 
with N q ~ q 2 at small g's but with the sign of the Fermi 
surface at large momenta. As shown previously in the 
limit of infinite Coulomb repulsioniS or in the low-density 
regime^ this WF has low-energy properties similar to 
one-dimensional Luttinger liquid conductors, where the 
absence of a jump in the momentum distribution is re- 
placed by a weaker singularity, yielding to 2kp and Akp 
power-law density correlations. It is important to em- 
phasize that, in the quantum case, the power-law corre- 
lations come from the large momentum singularity, that 
are absent in the classical CGMti4 Indeed, in the quan- 
tum state, the subleading corrections in the classical po- 
tential of Eq. 10 are very important and can actually 
turn the CGM exponential correlations to power laws in 
the plasma phase, and vice-versa in the confined phase. 
On the other hand, in the confined phase the Gaussian 
approximation is not adequate both at small and large 
momenta, see Fig- El Indeed, at small g's, the coefficient 
of the quadratic term is not simply given by the Gaus- 
sian approximation and, more importantly, the strong 
Jastrow factor washes out completely the singularities of 
iV° , leading to a smooth charge-structure factor, a gen- 
uine fingerprint of an insulating phase. 

In order to further characterize the two phases, we con- 
sider the quasiparticle weight 



\<$!N-l\ck, a \*N)? 

(*jv|*Jv>(*iv-i|*Ar-i>' 



(17) 



where |*jv) an d |*jv— l) are the WF with N and (N - 1) 
particles, Ck l(7 is the destruction operator of a particle of 
momentum fc and spin a. In particular, the wave function 
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FIG. 7: (Color online) Optimized Jastrow potential v q , mul- 
tiplied by \q\ 2 , for the Hubbard model as a function of \q\ [in 
the (1,1) direction] for different sizes of the cluster and ra- 



tios U/t. The arrow indicates 7r/T c 



the expected value 



at the classical transition point. 



with N — 1 particles is constructed from \^n)'- 



\^N-i) = Jc k>cr \V), 



(18) 



In a previous workj2& it was argued that the singular Jas- 
trow factor can induce non-Fermi liquid properties, and 
in particular a vanishing at the Fermi surface. In 
Fig. we report Z k for k — (tt/2,tt/2) and for differ- 
ent Jastrow strengths (3. We find that the quasiparticle 
weight vanishes with a power-law behavior: 



(19) 



both in the confined and in the plasma phase, with an 
exponent 9 that depends upon (3 and the type of the 
uncorrelated state. In the plasma phase, 9 varies contin- 
uously with the Jastrow strength and there is no ap- 
preciable dependence on the uncorrelated determinant. 
On the other hand, in the confined phase, the exponent 
is constant, i.e., 9 ~ 1/2 for the BCS state and 9 ~ 3/4 
for the free-electron state, and equal to the value found 
for the fully projected WF, as shown in Fig. It must 
be mentioned that, for the BCS state, 9 does not depend 
upon the value of the superconducting order parameter 
A (see Fig.0J, indicating the universal properties of the 
confined phase. 

Our results show that it is possible to describe a con- 
tinuous MIT in 2D electronic systems with a Jastrow 
correlated WF. We characterized both the metallic re- 
gion, with a zero-quasiparticle weight = 0, and the 
insulator, which can be continuously connected to the 
totally projected Gutzwiller WF. Of course, other sce- 
narios are possible for the MIT, e.g., the one proposed 
in the infinite-dimensional limit. Indeed, whenever the 
metallic phase has N q ~ c\q\ + dq 2 with c > 0, the MIT 
is not described by the functional form (|13l) of the Jas- 
trow: By approaching the transition from the metallic 
phase, we enter directly into the confined phase with a 
quadratic charge structure factor at small momenta, i.e., 
c — > with a large finite d at the critical point. In this 



case, in the metallic region, a less singular Jastrow factor 
v q ~ l/\q\ is expected, leading to a finite quasiparticle 
weight. 

Finally, we would like to comment about the possibility 
to stabilize the "Coulomb metal" phase in a microscopic 
model. For simplicity, let us consider the one-band Hub- 
bard model on the square lattice with nearest-neighbor 
hopping: 



H 



h.c. 



(20) 



where c\ a creates an electron with spin a at the site 

i, Tii ><7 = c\ a Ci, a is the density operator at the site i. 
In the following, we will consider a paramagnetic state, 
by taking a projected Fermi gas wave function, as de- 
scribed above, and by minimizing the variational energy 
for the determination of the Jastrow factor J . In partic- 
ular, by using the method described in Ref. l2ll we are 
able to optimize all the independent Jastrow parameters 
in the real space Vij (i.e., the Fourier transform of v q ). 
Within this approach, that neglects magnetic phases, we 
obtain a MIT for U c /t — 8.5 ± 0.5. In the weak cou- 
pling regime, for U < U c , we obtain a Fermi liquid with 
a finite quasiparticle weight, whereas at strong couplings, 
i.e., for U > U c , we have an insulating phase with a van- 
ishing Zk, see Fig. EI 22 Moreover, the calculation of the 
double occupancy D clearly indicates that the transition 
is continuous and the insulating phase still possess fi- 
nite charge fluctuations, see inset of Fig. EI As discussed 
above, in the metallic region, we find that v q ~ 1/|<7| (see 
Fig. UJ . Unfortunately, as soon as we enter in the insu- 
lating phase, lim g ^o v q \q\ 2 defines an effective (3 which 
is larger than the critical value for the KT transition, 
and, therefore, no evidence for the "Coulomb metal" is 
founds Indeed, we expect that the optimized Jastrow 
factor v q , containing subleading corrections with respect 
to Eq. I|13|l , will define a critical (3 very close to the value 
of the classical CGM, i.e., /3 C = 1/T C CGM = 4. There- 
fore, in light of the results of Fig. [7J the stabilization 
of the "Coulomb metal" seems to be very unlikely: Al- 
though there are large size effects around U c , we have a 
clear evidence that lim^o^M 2 — * for U < U c and 
lim 9 ^ v q \q\ 2 > 47T for U > U c (see Fig.Q. 

On the other hand, we can safely predict the occur- 
rence of the novel KT-like scenario described above in 
2D systems with long-range (logarithmic) interaction. In 
this case, the application of the Gaussian approximation 
for small interaction and our ansatz for the insulating 
phase imply the presence of a transition of the type con- 
sidered here. It is remarkable that the proposed pic- 
ture crucially depends on the long-range nature of the 
Coulomb interaction, recalling Mott's original idea. In 
this regard, it should be mentioned that the original Mott 
argument for a discontinuous metal-insulator transition, 
driven by the long-range Coulomb interaction, cannot 
be applied in 2D. In such case there always exists a 
bound state for two opposite charges interacting with the 
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screened Coulomb potential, so that, according to this ar- 
gument, no metallic phase with unbound charges is possi- 
ble. However, this is clearly an artifact of such mean-field 
argument, since for small interaction the random- phase 
approximation leads to an anomalous metallic state (see 
Ref.H3). 
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